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Abstract 


This article gives a brief history of the analysis and computation of the mathemat- 
ical constant n = 3.14159 ... , including a number of the formulas that have been 
used to compute n through the ages. Recent developments in this area are then 
discussed in some detail, including the recent computation of n to over six billion 
decimal digits using high-order convergent algorithms, and a newly discovered 
scheme that permits arbitrary individual hexadecimal digits of 7t to be computed. 
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Introduction 

The fascinating history of the constant we now know as 7r spans several millennia, 
almost from the beginning of recorded history up to the present day. In many ways this 
history parallels the advancement of science and technology in general, and of mathematics 
and computer technology in particular. An overview of this history is presented here in 
sections one and two. Some exciting recent developments are discussed in sections three 
and four. Section five explores the question of why this topic has such enduring interest. 
For further details of the history of tt up to about 1970, the reader is referred to Petr 
Beckmann’s readable and entertaining book [3]. A listing of milestones in the history of 
the computation of 7r is given in Tables 1 and 2. 

1. The Ancients 

In one of the earliest accounts (about 2000 BC) of 7 r, the Babylonians used the approx- 
imation 3| = 3.125. At this same time or earlier, according to an account in an ancient 
Egyptian document, Egyptians were assuming that a circle with diameter nine has the 
same area as a square of side eight, which implies n = ~p ~ 3.1604 . . .. Others of antiquity 
were content to use the simple approximation 3, as evidenced by the following passage 
from the Old Testament: 


Also, he made a molten sea of ten cubits from brim to brim, round in compass, 
and five cubits the height thereof; and a line of thirty cubits did compass it 
round about (I Kings 7:23; see also 2 Chron. 4:2). 

The first rigorous mathematical calculation of the value of 7r was due to Archimedes 
of Syracuse (ca. 250 BC), who used a geometrical scheme based on inscribed and circum- 
scribed polygons to obtain the bounds 3“ < 7r < 3^, or in other words 3.1408. . . < 7T < 
3.1428. .. [11]. No one was able to improve on Archimedes’ method for many centuries, 
although a number of persons used this general method to obtain more accurate approxi- 
mations. For example, the astronomer Ptolemy, who lived in Alexandria in 150 AD, used 
the value 3^ = 3.141666..., and the fifth century Chinese mathematician Tsu Chung- 
Chih used a variation of Archimedes’ method to compute 7 r correct to seven digits, a level 
not obtained in Europe until the 1500s. 


2. The Age of Newton 

As in other fields of science and mathematics, little progress was made in the quest for 
7 r during the dark and middle ages, at least in Europe. The situation was somewhat better 
in the East, where Al-Kashi of Samarkand computed 7r to 14 places about 1430. But in 
the 1600s, with the discovery of calculus by Newton and Leibniz, a number of substantially 
new formulas for it were discovered. One of them can be easily derived by recalling that 
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Substituting x = 1 gives the well-known Gregory-Leibniz formula 


7r / 4 = 1 -1/3 + 1/5-1/7 + 1/9-1/11 +... 


Regrettably, this series converges so slowly that hundreds of terms would be required to 
compute the numerical value of 7 r to even two digits accuracy. However, by employing the 
trigonometric identity 

7r/4 = tan _1 (l/2) + tan _1 (l/3) 

(which follows from the addition formula for the tangent function), one obtains 



which converges much more rapidly. An even faster formula, due to Machin, can be 
obtained by employing the identity 

?r/4 = 4tan _1 (l/5) — tan~ 1 (l/239) 


in a similar way. Shanks used this scheme to compute n to 707 decimal digits accuracy in 
1873. Alas, it was later found that this computation was in error after the 527-th decimal 
place. 

Newton discovered a similar series for the arcsine function: 


sin 1 x 


1 • x 3 1 • 3 * x 5 1 ■ 3 • 5 ■ x 7 

= * + — r+ ; z + — — : — - — r + • 


2*3 


2-4*5 2 • 4 • 6 • 7 


7r can be computed from this formula by noting that 7r /6 = sin 1 (l/2). An even faster 
formula of this type is 


2^1 + 24 M ! I L 

4 \3 • 2 3 5 - 2 s 7 ■ 2 7 9 • 2 9 


Newton himself used this particular formula to compute 7 r. He published 15 digits, but 
later sheepishly admitted, “I am ashamed to tell you how many figures I carried these 
computations, having no other business at the time.” 

In the 1700s the mathematician Euler, arguably the most prolific mathematician in 
history, discovered a number of new formulas for 7 r. Among these are 


7i 1111 

7F ~ 1 + 22 + 32 + 42 + 52 + 

7 r 4 _ 1111 , 

90 _ 1 + ¥ + ¥ + 4 ^ + 5 1 + ' " 


A related, more rapidly convergent series is 


6 


OO 


3£- 

m = 1 m 



3 



These formulas aren’t very efficient for computing 7 r, but they have important theoreti- 
cal implications and have been the springboard for notable research questions, such as the 
Riemann zeta function hypothesis, that continue to be investigated to this day. 

One motivation for computations of i r during this time was to see if the decimal expan- 
sion of 7r repeats, thus disclosing that i r is the ratio of two integers (although hardly anyone 
in modern times seriously believed that it was rational). This question was conclusively 
settled in the late 1700s, when Lambert and Legendre proved that 7r is irrational. Some 
still wondered whether 7 r might be the root of some algebraic equation with integer coeffi- 
cients (although as before few really believed that it was). This question was finally settled 
in 1882 when Lindemann proved that 7r is transcendental. Lindemann’s proof also settled 
once and for all, in the negative, the ancient Greek question of whether the circle could 
be squared with ruler and compass. This is because constructible numbers are necessarily 
algebraic. 

In the annals of 7T, the nineteenth century came to a close on an utterly shameful note. 
Three years prior to the turn of the century, one Edwin J. Goodman, M.D. introduced into 
the Indiana House of Representatives a bill that would introduce “new Mathematical truth” 
and enrich the state, which would profit from the royalties ensuing from this discovery. 
Section two of the bill included the passage 

“disclosing the fourth important fact that the ratio of the diameter and cir- 
cumference is as five-fourths to four;” 

Thus one of Goodman’s new mathematical “truths” is that 7r = y = 3.2. In spite of this 
and numerous other absurd statements, the Indiana House passed the bill unanimously 
on Feb. 5, 1897. The bill then passed a Senate committee, and would have been enacted 
into law had it not been for the last-minute intervention of Prof. C. A. Waldo of Purdue 
University, who happened to hear some of the deliberation while on other business. 

3. The Twentieth Century 

With the development of computer technology in the 1950s, 7 r was computed to thou- 
sands and then millions of digits, in both decimal and binary bases (see for example [17]). 
These computations were facilitated by the discovery of some advanced algorithms for per- 
forming the required high-precision arithmetic operations on a computer. For example, in 
1965 it was found that the newly-discovered fast Fourier transform (FFT) could be used 
to perform high-precision multiplications much more rapidly than conventional schemes. 
These methods dramatically lowered the computer time required for computing 7r and 
other mathematical constants to high precision. See [1], [7] and [8] for a discussion of some 
of these techniques. 

In spite of these advances, until the 1970s all computer evaluations of 7r still employed 
classical formulas, usually a variation of Machin’s formula. Some new infinite series formu- 
las were discovered by the Indian mathematician Ramanujan around 1910, but these were 
not well known until quite recently when his writings were widely published. One of these 
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is the remarkable formula 


1 

7T 


2y/2 ~ 

9801 to 


(4Jb)!(1103 + 26390*) 
(fc!) 4 396 4 * 


Each term of this series produces an additional eight correct digits in the result. Gosper 
used this formula to compute 17 million digits of n in 1985. 

While Ramanujan’s series is considerably more efficient than the classical formulas, 
it shares with them the property that the number of terms one must compute increases 
linearly with the number of digits desired in the result. In other words, if one wishes to 
compute 7 r to twice as many digits, then one must evaluate twice as many terms of the 
series. 

In 1976 Eugene Salamin [16] and Richard Brent [8] independently discovered a new 
algorithm for 7r, which is based on the arithmetic-geometric mean and some ideas originally 
due to Gauss in the 1800s (although for some reason Gauss never saw the connection to 
computing 7 r). This algorithm produces approximations that converge to tt much more 
rapidly than any classical formula. The Salamin-Brent algorithm may be stated as follows. 
Set a 0 = 1, b 0 = l/y/2 and s 0 = 1/2. For k = 1,2, 3, • • • compute 


CLk 


h 

Ck 

•sjt 

Pk 


a k~ 1 + frfc-1 


\J Gtjfc-l6fc-l 

2 l2 
a k ~ °k 

$k — 1 2 Ck 

2 a| 

Sk 


Then pk converges quadratically to 7r. This means that each iteration of this algorithm 
approximately doubles the number of correct digits. To be specific, successive iterations 
produce 1, 4, 9, 20, 42, 85, 173, 347 and 697 correct digits of n. Twenty-five iterations are 
sufficient to compute 7 r to over 45 million decimal digit accuracy. However, each of these 
iterations must be performed using a level of numeric precision that is at least as high as 
that desired for the final result. 

The Salamin-Brent algorithm requires the extraction of square roots to high precision, 
operations not required, for example, in Machin’s formula. High-precision square roots 
can be efficiently computed by means of a Newton iteration scheme that employs only 
multiplications, plus some other operations of minor cost, using a level of numeric precision 
that doubles with each iteration. The total cost of computing a square root in this manner 
is only about three times the cost of performing a single full-precision multiplication. 
Thus algorithms such as the Salamin-Brent scheme can be implemented very rapidly on a 
computer. 

Beginning in 1985, two of the present authors (Jonathan and Peter Borwein) discovered 
some additional algorithms of this type [5, 6, 7]. One is as follows. Set a 0 = 1/3 and 
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so = (v^ - l)/2. Iterate 

3 

1+2(1 -.,^1/3 

n + 1 - l 

2 

r k+l a k ~~ 3 ( t + + 1 _ 

Then 1 /a* converges cubically to 7r — each iteration approximately triples the number of 
correct digits. A quartic algorithm is as follows: Set a 0 = 6 — 4\/2 and j/ 0 = \/2 — 1. Iterate 

i + (i - ytY 14 

Ofc(l + Vk+ 1 ) 4 - 2 2k+3 y k+1 (l + Vk+ 1 + y* +1 ) 

Then ajt converges quartically to 1/tt. This particular algorithm, together with the Salamin- 
Brent scheme, has been employed by Yasumasa Kanada of the University of Tokyo in 
several computations of n over the past ten years or so. In the latest of these computations, 
Kanada computed over 6.4 billion decimal digits on a Hitachi supercomputer. This is 
presently the world’s record in this arena. 

More recently it has been further shown that there are algorithms that generate m-th 
order convergent approximations to 7 r for any m. An example of a nonic (ninth-order) 
algorithm is the following: Set a 0 = 1/3, r 0 = (\/3 - l)/2, s 0 = (1 - r^) 1 / 3 . Iterate 

1 + 2 r k 

[9r*(l +r k + r'l)Y /3 
t 2 + tu + u 2 
27(1 + s k + s 2 .) 

v 

ma k + 3 2fc_1 (l — m) 

(1-nO 3 

( t + 2u)v 

( i-.D' /3 

Then 1 ja k converges nonically to 7 r. It should be noted however that these higher order 
algorithms do not appear to be faster as computational schemes than, say the Salamin- 
Brent or the Borwein quartic algorithms. In other words, although fewer iterations are 
required to achieve a given level of precision in the higher-order schemes, each iteration is 
more expensive. 

A comparison of actual computer run times for various ir algorithms is shown in Figure 
1. These run times are for computing 7 r in binary to various precision levels on an IBM 
RS6000/590 workstation. The abscissa of this plot is in hexadecimal digits — multiply 
these numbers by four to obtain equivalent binary digits, or by log 10 (16) = 1.20412... 


t = 
u = 
v = 

m = 

a k + 1 = 

•Sfc + l = 
r *+i = 


Vk+i = 
«k+i = 


n t+i = 

■Sfc+l = 
a k+\ = 
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Run time for computing various n algorithms 



Time in seconds 

3000 




Number of hex digits 


Gregory-Leibniz 

Machin 

Brent-Salamin no FFT 
Brent-Salamin FFT 
Borwein quartic FFT 
Borwein cubic FFT 
Borwein Garvan nonic FFT 


Figure 1: Run times for Computing Pi 


to obtain equivalent decimal digits. Other implementations on other systems may give 
somewhat different results — for example, in Kanada’s recent computation of i r to over six 
billion digits, the quartic algorithm ran somewhat faster than the Salamin- Brent algorithm 
(116 hours versus 131 hours). But the overall picture from such comparisons is unmistak- 
able: the modern schemes run many times faster than the classical schemes, especially 
when implemented using FFT-based arithmetic. 

David and Gregory Chudnovsky of Columbia University have also done some very- 
high precision computations of 7r in recent years, alternating with Kanada for the world’s 
record. Their most recent computation (1994) produced over four billion digits of 7 r [9]. 
They did not employ a high-order convergent algorithm, such as the Salamin-Brent or 
Borwein algorithms, but instead utilized the following infinite series (which is in the spirit 
of Ramanujan’s series above): 

1 ~ ( — 1)* (6^)1(13591409 + 545140134*) 

7 ” 12 (3£-)!(fc!) 3 640320 3 *+ 3 / 2 


Each term of this series produces an additional 14 correct digits. The Chudnovskys imple- 
mented this formula with a very clever scheme that enabled them to utilize the results of a 
certain level of precision to extend the calculation to even higher precision. Their program 
was run on a home-brew supercomputer that they have assembled using private funds. An 


7 




interesting personal glimpse of the Chudnovsky brothers is given in [14]. 

4. Computing Individual Digits of n 

At several junctures in the history of 7 r, it was widely believed that virtually everything 
of interest with regards to this constant had been discovered, and in particular that no 
fundamentally new formulas for 7 r lay undiscovered. This sentiment was even suggested in 
the closing chapters of Beckmann’s 1971 book on the history of 7 r [3, pg. 172]. Ironically, 
the Salamin- Brent algorithm was discovered only five years later. 

A more recent reminder that we have not come to the end of humanity’s quest for 
knowledge about 7r came with the discovery of the Rabinowitz- Wagon “spigot” algorithm 
for 7r in 1990 [15]. In this scheme, successive digits of 7r (in any desired base) can be 
computed with a relatively simple recursive algorithm based on the previously generated 
digits. Multiple precision computation software is not required, so that this scheme can be 
easily implemented on a personal computer. 

Note however that this algorithm, like all of the other schemes mentioned above, still has 
the property that in order to compute the d-th digit of 7r, one must first (or simultaneously) 
compute each of the preceding digits. In other words, there is no “shortcut” to computing 
the d- th digit with these formulas. Indeed, it has been widely assumed in the field (although 
never rigorously proven) that the computational complexity of computing the d - th digit is 
not significantly less than that of computing all of the digits up to and including the d - th 
digit. This may still be true, although it is probably very hard to prove. Another common 
feature of the previously known 7T algorithms is that they all appear to require substantial 
amounts of computer memory, amounts that typically grow linearly with the number of 
digits generated. 

Thus it was with no small surprise that a novel scheme was recently discovered for 
computing individual hexadecimal digits of 7r [2]. In particular, this algorithm (1) pro- 
duces the d - th hexadecimal (base 16) digit of 7r directly, without the need of computing 
any previous digits; (2) is quite simple to implement on a computer; (3) does not require 
multiple precision arithmetic software; (4) requires very little memory; and (5) has a com- 
putational cost that grows only slightly faster than the index d . For example, the one 
millionth hexadecimal digit 7r can be computed in only a minute or two on a current RISC 
workstation or high-end personal computer. This algorithm is not fundamentally faster 
than other known schemes for computing all digits up to some position d, but its elegance 
and simplicity are nonetheless of considerable interest. 

This scheme is based on the following remarkable new formula for i r: 

^ 1 / 4 2 1 1 \ 

^ ~ “ 16 1 V 82 -F 1 8 i + 4 8?: + 5 8i + 6J 

The proof of this formula is not very difficult. First note that for any k < 8, 
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Babylonians 

2000? BCE 

1 

3.125 (3f) 

Egyptians 

2000? BCE 

1 

3.16045 (4(§) 2 ) 

China 

1200? BCE 

1 

3 

Bible ( 1 Kings 7:23) 

550? BCE 

1 

3 

Archimedes 

250? BCE 

3 

3.1418 (ave.) 

Hon Han Shu 

130 AD 

1 

3.1622 (= V^0 ?) 

Ptolemy 

150 

3 

3.14166 

Chung Hing 

250? 

1 

3.16227 (vTO) 

Wang Fau 

250? 

1 

3.15555 (ff ) 

Liu Hui 

263 

5 

3.14159 

Siddhanta 

380 

3 

3.1416 

Tsu Ch’ung Chi 

480? 

7 

3.1415926 

Aryabhata 

499 

4 

3.14156 

Brahmagupta 

640? 

1 

3.162277 (= v/T0) 

Al-Khowarizmi 

800 

4 

3.1416 

Fibonacci 

1220 

3 

3.141818 

Al-Kashi 

1429 

14 


Otho 

1573 

6 

3.1415929 

Viete 

1593 

9 

3.1415926536 (ave.) 

Romanus 

1593 

15 


Van Ceulen 

1596 

20 


Van Ceulen 

1615 

35 


Newton 

1665 

16 


Sharp 

1699 

71 


Seki 

1700? 

10 


Kamata 

1730? 

25 


Machin 

1706 

100 


De Lagny 

1719 

127 

(112 correct) 

Takebe 

1723 

41 


Matsunaga 

1739 

50 


Vega 

1794 

140 


Rutherford 

1824 

208 

(152 correct) 

Strassnitzky and Dase 

1844 

200 


Clausen 

1847 

248 


Lehmann 

1853 

261 


Rutherford 

1853 

440 


Shanks 

1874 

707 

(527 correct) 


Table 1: History of x Calculations (Pre 20th Century) 
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Ferguson 

Ferguson 

Ferguson and Wrench 

Smith and Wrench 

Reitwiesner et al. (ENIAC) 

Nicholson and Jeenel 

Felton 

Genuys 

Felton 

Guilloud 

Shanks and Wrench 
Guilloud and Filliatre 
Guilloud and Dichampt 
Guilloud and Bouyer 
Miyoshi and Kanada 
Guilloud 
Tamura 

Tamura and Kanada 
Tamura and Kanada 
Kanada, Yoshino and Tamura 
Ushiro and Kanada 
Gosper 
Bailey 

Kanada and Tamura 

Kanada and Tamura 

Kanada, Tamura, Kuho, et. al 

Kanada and Tamura 

Chudnovskys 

Chudnovskys 

Kanada and Tamura 

Kanada and Tamura 

Chudnovskys 

Chudnovskys 

Chudnovskys 

Takahashi and Kanada 

Kanada 

Kanada 


1946 

620 

Jan. 1947 

710 

Sep. 1947 

808 

1949 

1,120 

1949 

2,037 

1954 

3,092 

1957 

7,480 

Jan. 1958 

10,000 

May 1958 

10,021 

1959 

16,167 

1961 

100,265 

1966 

250,000 

1967 

500,000 

1973 

1,001,250 

1981 

2,000,036 

1982 

2,000,050 

1982 

2,097,144 

1982 

4,194,288 

1982 

8,388,576 

1982 

16,777,206 

Oct. 1983 

10,013,395 

1985 

17,526,200 

Jan. 1986 

29,360,111 

Sep. 1986 

33,554,414 

Oct. 1986 

67,108,839 

Jan. 1987 

134,217,700 

Jan. 1988 

201,326,551 

May 1989 

480,000,000 

Jun. 1989 

525,229,270 

Jul. 1989 

536,870,898 

Nov. 1989 

1,073,741,799 

Aug. 1989 

1,011,196,691 

Aug. 1991 

2,260,000,000 

May 1994 

4,044,000,000 

Jun. 1995 

3,221,225,466 

Aug. 1995 

4,294,967,286 

Oct. 1995 

6,442,450,938 


Table 2: History of it Calculations (20th Century) 
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Thus we can write 


^ J_ /_4 2 1 J_\ 

16' \8i + 1 8* + 4 8t + 5 8> + 6/ 

4^/2 — 8x 3 — 4\/2x 4 — 8x 5 
Jo 1 — x 8 

which on substituting y := \/2x becomes 


Jo y 4 — 


2y 3 + 4y — 4 


J -r^ d y-J 0 

Jo y ^ — Z Jo 


4y - 8 
y 2 - 2j/ + 2 


reflecting a partial fraction decomposition of the integral on the left-hand side. 

However, this derivation is dishonest, in the sense that the actual route of discovery was 
much different. This formula was actually discovered not by formal reasoning, but instead 
by numerical searches on a computer using the “PSLQ” integer relation finding algorithm 
[10]. Only afterwards was a rigorous proof found. 

A similar formula for 7 r 2 (which also was first discovered using the PSLQ algorithm) is 
as follows: 


~ 1 [ 16 16 8 16 
^ I6 7 [( 87 + l ) 2 ~~ ( 8 * + 2) 2 ~ ( 8 * + 3) 2 ( 8 a + 4) 2 

4 4 2 

“ ( 8 * + 5) 2 “ ( 8 a + 6) 2 + ( 8 * + 7 ) 2 . 


Formulas of this type for a few other mathematical constants are given in [2]. 

Computing individual hexadecimal digits of it using the above formula crucially relies 
on what is known as the binary algorithm for exponentiation, wherein one evaluates x n 
by successive squaring and multiplication. This reduces the number of multiplications 
required to less than 21og 2 (n). According to Knuth, this technique dates back at least to 
200 B.C [13]. In our application, we need to obtain the exponentiation result modulo a 
positive integer c. This can be efficiently done with the following variant of the binary 
exponentiation algorithm, wherein the result of each multiplication is reduced modulo c: 
To compute r = b n mod c, first set t to be the largest power of two < n, and set r = 1. 
Then 

A: if n > t then r *— br mod c; n +— n — <; endif 
t t/2 

if t > 1 then r *— r 2 mod c; go to A; endif 

Here “mod” is used in the binary operator sense, namely as the binary function defined by 
x mod y := x - [xjy]y. Note that the above algorithm is entirely performed with positive 
integers that do not exceed c 2 in size. 

Consider now the first of the lour sums in the formula above for 7 r: 
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First observe that the hexadecimal digits of *S’i beginning at position d+ 1 can be obtained 
from the fractional part of I6^5‘i. Then we can write 


frac(16^S'i) 


OG 


E 


16 rf - fc 
8 k + 1 


mod 1 


d 


E 


\6 d - k mod 8k + 1 
8k + 1 


mod 1 


+ 


OO 


E 

k=d+\ 


l6 d-k 

8 k + 1 


mod 1 


For each term of the first summation, the binary exponentiation scheme can be used to 
rapidly evaluate the numerator. In a computer implementation this can be done using ei- 
ther integer or 64-bit floating-point arithmetic. Then floating-point arithmetic can be used 
to perform the division and add the quotient to the sum mod 1. The second summation, 
where the exponent of 16 is negative, may be evaluated as written using floating-point 
arithmetic. It is only necessary to compute a few terms of this second summation, just 
enough to insure that the remaining terms sum to less than the “epsilon” of the floating- 
point arithmetic being used. The final result, a fraction between 0 and 1, is then converted 
to base 16, yielding the ( d + l)-th hexadecimal digit, plus several additional digits. Full 
details of this scheme, including some numerical considerations, as well as analogous for- 
mulas for a number of other basic mathematical constants, can be found in [2]. Sample 
implementations of this scheme in both Fortran and C are available from the web site 
http : //www . cecm. sfu . ca/personal/pborwein/. 

As the reader can see, there is nothing very sophisticated about either this new formula 
for 7 r, its proof, or the scheme just described to compute hexadecimal digits of 7 r using it. 
In fact, this same scheme can be used to compute binary (or hexadecimal) digits of log(2) 
based on the formula 


l°g(2) = 


E 


*:=i 


i 

W' 


which has been known for centuries. Thus it is frankly astonishing that these methods 
have lain undiscovered all this time. There seems to be no fundamental reason that Euler, 
for example, could not have discovered them. The only advantage that today’s researchers 
have in this regard is advanced computer technology. Along this line, Table 3 gives some 
hexadecimal digits of 7 r computed using the above scheme. 

One question that immediately arises in the wake of this discovery is whether or not 
there is a formula of this type and an associated computational scheme to compute indi- 
vidual decimal digits of 7T. Alas, no decimal scheme for 7r is known at this time, although 
there is for certain constants such as log(9/10) — see [2]. On the other hand, there is 
not yet any proof that a decimal scheme for n cannot exist. This question is currently 
being actively pursued by researchers. Based on some numerical searches using the PSLQ 
algorithm, it appears that there are no simple formulas of the above form for it with 10 
in the place of 16. This of course does not rule out the possibility of completely different 
formulas that nonetheless permit rapid computation of individual decimal digits of 7 r. 
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Position 

Hex Digits Beginning 
At This Position 

10 6 

26C65E52CB4593 

10 7 

17AF5863EFED8D 

10 8 

ECB840E21926EC 

10 9 

85895585A0428B 

io 10 

921C73C6838FB2 


Table 3: Hexadecimal Digits of i r 


5. Why? 

Certainly there is no need for computing i r to millions or billions of digits in practical 
scientific or engineering work. A value of 7 r to 40 digits would be more than enough 
to compute the circumference of the Milky Way galaxy to an error less than the size of 
a proton. There are certain scientific calculations that require intermediate calculations 
to be performed to significantly higher precision than required for the final results, but 
it is doubtful than anyone will ever need more than a few hundred digits of 7r for such 
purposes. Values of 7r to a few thousand digits are sometimes employed in explorations of 
mathematical questions using a computer, but we not aware of any significant number of 
applications beyond this level. 

One motivation for computing digits of it is that these calculations are excellent tests 
of the integrity of computer hardware and software. This is because if even a single error 
occurs during a computation, almost certainly the final result will be in error. On the 
other hand, if two independent computations of digits of it agree, then most likely both 
computers performed billions or even trillions of operations flawlessly. For example, in 
1986, a 7 r-calculating program detected some obscure hardware problems in one of the 
original Cray-2 supercomputers [1]. 

The challenge of computing it has also stimulated research into advanced computational 
techniques. For example, some new techniques for efficiently computing linear convolutions 
and fast Fourier transforms (FFTs), which have applications in many areas of science and 
engineering, had their origins in efforts to accelerate computations of it. 

Beyond immediate practicality, decimal and binary expansions of it have long been 
of interest to mathematicians, who have still not been able to resolve the question of 
whether the expansion of it is normal [18]. In particular, it is widely suspected that the 
decimal expansions of it, e, \/2, v^lO, and many other mathematical constants all have the 
property that the limiting frequency of any digit is one tenth, and the limiting frequency 
of any n-long string of decimal digits is 10 _ ” (and similarly for binary expansions). Such a 
guaranteed property could, for instance, be the basis of a reliable pseudo-random number 
generator for scientific calculations. Unfortunately, this assertion has not been proven in 
even one instance. Thus there is a continuing interest in performing statistical analyses on 
the expansions of these numbers to see if there is any irregularity that would suggest this 
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assertion is false. So far, such studies of high-precision values of n have not disclosed any 
irregularities. Along this line, new formulas and schemes for computing digits of 7 r, such 
as the one described in section four, are of interest because some of these may suggest new 
approaches to answering the normality question. 

Finally, there is a more fundamental motivation for computing 7 r, which should be 
familiar to anyone who has scaled a lofty mountain or competed in a major sporting event: 
“it is there ’ 1 — it is easily the most famous of the basic constants of mathematics. Thus 
as long as there are humans (and computers) we will doubtless have ever-more impressive 
computations of 7 r. 

Conclusion 

The constant 7 r has repeatedly surprised humanity with new and often unanticipated 
results. If anything, the discoveries of this century have been even more startling, with 
respect to the previous state of knowledge, than those of past centuries. Thus we conclude 
that even more surprises lurk in the depths of undiscovered knowledge regarding this famous 
constant. We thus look forward to what the future has to bring. 
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